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We report the first measurement of the diffusion coefficients of francium and rubidium ions im- 
planted in a yttrium foil. We developed a methodology, based on laser spectroscopy, which can be 
applied to radioactive and stable species, and allows us to directly take record of the diffusion time. 
Francium isotopes are produced via fusion-evaporation nuclear reaction of a 100 MeV 18 O beam 
on a Au target at the Tandem XTU accelerator facility in Legnaro, Italy. Francium is ionized at 
the gold-vacuum interface and Fr + ions are then transported with a 3 keV electrostatic beamline 
to a cell for neutralization and capture in a magneto-optical trap (MOT). A Rb + beam is also 
available, which follows the same path as Fr + ions. The accelerated ions are focused and implanted 
in a 25 /im thick yttrium foil for neutralization: after diffusion to the surface, they are released as 
neutrals, since the Y work function is lower than the alkali ionization energies. The time evolution 
of the MOT and the vapor fluorescence signals are used to determine diffusion times of Fr and Rb 
in Y as a function of temperature. 



I. INTRODUCTION 

In many experiments with radioactive atoms, elements 
are produced as ions and then neutralized. The most im- 
portant requirement for a neutralization system is the 
fast release of neutral atoms, with respect to their ra- 
dioactive decay time. In particular, in our francium ex- 
periment at INFN's National Laboratories in Legnaro, 
Italy [l|, we produce Fr ions that have to be neutralized 
before accumulation in a magneto-optical trap. The sys- 
tem must be very efficient, as a large sample of Fr atoms 
is requested for nuclear decay, atomic parity violation 
and permanent electric dipole studies. 

Our apparatus follows the scheme of the Stony Brook 
experiment @, where many measurements on Fr have 
been performed so far. Francium isotopes are produced 
at the Tandem accelerator facility in Legnaro via the 
fusion-evaporation reaction 

197 Au+ 18 0— » 21B - x ft + sn. 

With an energy of the primary beam around 100 MeV, 
we are able to produce Fr isotopes in the mass number 
range 208-211. These isotopes have a lifetime that is long 
enough for laser cooling and trapping. The system is at 
the moment optimized for maximum production of 210 Fr, 
which has a half-life of 191 s. Francium atoms are ionized 
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at the surface of the gold target and transported through 
an electrostatic beamline at an energy of 3 keV towards 
a spectroscopic cell designed for laser cooling. The ionic 
beam is focused on a 99.96% pure yttrium foil, placed 
inside the cell on the opposite side with respect to the 
entrance aperture. In order to be able to detect the trap, 
we have to optimize all the processes, from production 
to laser trapping. With only a few days of beamtime 
per year, it is not convenient to test the setup with fran- 
cium: we decided to use stable rubidium, always available 
in high quantities. To this purpose a Rb dispenser was 
placed in the reaction chamber near the target: Rb atoms 
that arrive on the gold surface are ionized and injected 
in the electrostatic line, following the same path as Fr. 
A sketch of the apparatus is shown in Fig. [TJ 

Yttrium has a work function <f> = 3.1 eV, that is lower 
than Fr and Rb ionization potentials, respectively 2p r = 
4.08 eV and Ijy-, = 4.18 eV. The probability of release in 
neutral form at the Y surface is close to unity, according 
to the Saha-Langmuir equation 

n+ g+ (4>- 1\ 
— = — ex P T^F > 

where n + and n a are the number of desorbed ions and 
neutral atoms, g + and g a are the statistical weights due 
to degeneracy of the ion and atom ground states respec- 
tively (jt- = I for alkalis), and ks is Boltzmann's con- 
stant. This ratio is - 6 x 10~ 6 at 1000 K. 

In spite of the fact that high temperature enhances dif- 
fusion, we fixed an upper limit of 1050 K in order to pre- 
serve the cell coating [3[ and hence not to compromise the 
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FIG. 1: Scheme of the experimental setup (the two anti- 
Helmholtz coils for magneto optical trapping are not shown) . 
In the inset, an image of the cold Fr trap is reported. The 
image is not to scale. 



trapping efficiency. This is a limitation because melting 
point of Y is T m =1799 K, and because the total release 
fraction grows with the ratio T/T m 0]. In this context, 
our goal was to check whether the efficiency of the neu- 
tralizer was high enough at a chosen temperature. So far 
the release efficiency of radioactive atoms embedded in 
the neutralizcr has been measured by monitoring their 
nuclear decay [3]. We used atomic laser spectroscopy to 
observe the time evolution of the released neutrals and 
deduce the diffusion time. We also took advantage of the 
extraordinary sensitivity given by our magneto-optical 
trap, which allows us to observe as few as 50 atoms. 



II. 1-D DIFFUSION MODEL 

If an ion is implanted at a given depth inside the Y foil, 
then it has to diffuse towards the surface to be released 
in atomic form. In order to achieve a good release effi- 
ciency, the time necessary to complete this process has 
to be small compared with the radioactive lifetime of the 
investigated atom. Since the transverse dimension of the 
ionic beam is much larger than the penetration depth, 
we can use the one-dimensional diffusion equation for the 
concentration N(x, t) of the diffusing species in the Y foil 
at the distance x from the surface: 



8N(x,t) 
di 



D 



d 2 N{x,t) 



T r N(x,t)+cp(x,t). 



r r is the radioactive decay rate (T r = for Rb), D the 
diffusion coefficient and tp(x,t) is the incoming current 
of ions impinging on the Y foil. Note that ip{x, t) is not 
directly observable: it has to be deduced from the mea- 
surement of its integral I(t) (the total current) on the 
neutralizer volume and from the presumed implantation 
distribution, discussed in Sections III Al and III Bl 

Empirically, we expect that the dependence of D on 



the temperature follows the Arrhenius' law 



D = Doc exp 



knT 



E a is the activation energy and is the theoretical 
asymptotic value of D for large temperatures. The char- 
acteristic diffusion time is given by 



Td 



AD ' 



where d is the mean implantation depth of ions in the Y 
foil. 



A. Pulsed regime 

In the case of No ions all implanted at the same dis- 
tance d from the surface at the time t = 0, i.e. N (x, 0) = 
NqS(x — d, 0), the diffusion equation with <p(x, t) = has 
an analytical solution: the flux of neutral atoms released 
in the cell at the time t is given by Q 



Fo(t) 




exp (-yj exp (-IV). 



(1) 



For an implantation distribution with a finite extent, the 
function can be different, specially for t < t^; in a rough 
approximation, this difference can be taken into account 
with a correction factor, which is expected to be of the 
order of unity in the case that r r <C t^ 1 Q. For instance, 
Melconian et al. take the case of a distribution modelled 
by a Gaussian times a linear term, with a characteris- 
tic implantation depth d Q. In this case, the function 
becomes 



Fo{t) 
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1 



No 

2 T d (l + f/ Td )3/2 



exp(— T r t). 



(2) 



B. Continuous regime 

If a constant total ionic current Iq is sent to the neu- 
tralizer starting at the time t = 0, for stable isotopes 
(r r = 0) we see that the neutralized released current 
takes the form 



F(t) = 



1 T m 



t') 



No 



dt' = Io 




(3) 



In the case of the implantation distribution considered 
by Melconian et at, we have 



F(t) = Io 



1 



VI + t/T d 



(4) 



We will see in the following that we can operate in two 
modes. If we let ions accumulate in cold yttrium for a 
certain time and then we turn on the neutralizer (pulsed 
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mode), all the implanted ions are released according to 
Eqs. © or ©. Instead, if we turn on the ionic beam 
with the neutralizer already working (continuous mode) , 
we are in the case described by Eqs. ([3]) or ©. 

In our model of the release process, we neglect the 
contribution of desorption to the release time. We can 
reasonably assume that, due also to the accuracy of such 
a model, discussed in Section [IV Bl desorption contribu- 
tion is well included in our uncertainty. In any case, as 
also discussed in Ref. [1], the diffusion coefficients de- 
rived from our data analysis can be properly considered 
as upper limits. 



In practice, provided that / varies with a time scale much 
longer than 1/C, we obtain that the detected number of 
trapped atoms and the total number of atoms in the va- 
por phase are proportional to the released current at any 
time. Therefore, the measurement of the trap signal or of 
the fluorescence signal directly gives all the information 
about diffusion inside the neutralizer. This constitutes 
an innovative way of measuring the diffusion coefficient, 
with respect to other experiments that measure the re- 
lease efficiency by detecting the radioactivity from the 
residual ions in the neutralizer. 



III. MAGNETO-OPTICAL TRAPPING 
DYNAMICS 

After being released, atoms can be either directly de- 
tected in the vapor phase (rubidium), or they can be 
first collected in a small volume (about 1 mm 3 ) in or- 
der to enhance our sensitivity (rubidium and francium). 
The confinement is provided by a magneto-optical trap in 
the standard configuration: six trapping and repumping 
laser beams respectively tuned on the D2 and Dl lines, 
and a constant gradient magnetic field provided by two 
anti-Helmholtz coils. 

The total number of atoms in the cell, when the six 
laser beams and the magnetic field are present, is the 
sum of two components: the number of trapped atoms 
Nt and the number of atoms in the vapor phase N v which 
is supposed to be at thermal equilibrium with the cell 
walls. The time evolution of iVj_and N v is well described 
by the coupled rate equations 



N t = LN V - CN t - T r iV t 



-(£ + W)N V + CN t - T r N v + f , 



(5) 



where L is the rate related to the loading process (the 
product LN V represents the number of atoms per unit 
time loaded from the vapor to the MOT), C is the rate de- 
scribing the collisional loss of atoms from the MOT (CNt 
represents the number of trapped atoms per unit time 
which go back from trap to vapor phase due to collisions 
with background gas in the cell) and W is the rate char- 
acterizing the loss of atoms from the vapor phase (WN V 
represents the number of atoms per unit time definitively 
lost to the vacuum system or by chemisorption to the cell 
walls). / is the current of atoms released by the neu- 
tralizer and suitable for trapping. We did not report the 
quadratic terms due to collisions between trapped atoms, 
which are negligible. We can find the quasi-stationary 
state solution Nt — N v ~ of Eqs.© under the condi- 
tion /// -C C,W (slowly varying incoming current hy- 
pothesis): 



Nt(t) 



fit) 



(w + r r )(c + r r ) + Lr r 

C + T r (6) 

Nv{t) = (w + r r )(c + r r ) + Lr r /(<) 



IV. MEASUREMENTS METHODS 
A. Vapor fluorescence and MOT detection 

At room temperature the Rb vapor density is high 
enough to be directly detected through resonant exci- 
tation and looking at the atomic fluorescence. In order 
to minimize changes in the apparatus we made measure- 
ments by only switching off the trapping magnetic field, 
but keeping in place the six laser beams. The laser fre- 
quency is scanned across the whole profile of the D2 Rb 
line to search for the maximum signal. Fluorescence light 
is in this case collected on a Si photodiode. We tried dif- 
ferent spectroscopic schemes and we obtained the best 
signal-to-noise ratio with the one described above. In 
the case of francium, it is not possible to directly observe 
the atomic vapor fluorescence, because the signals are ex- 
tremely weak. The solution is given by the MOT, that 
allows to accumulate Fr (and also Rb) atoms in a small 
volume at low temperature. In fact the signal increases 
by about two orders of magnitude because of Doppler 
broadening elimination and, at low densities, far from 
saturated vapor pressures, by several orders of magni- 
tude because of higher densities. As already stressed, 
we are able to detect traps of only a few tens of atoms. 
Atomic fluorescence from trapped atoms is detected by a 
CCD cooled camera (Fig.©. A dedicated software per- 
forms background subtraction, to reduce the contribution 
of scattered light to the noise, and real time correction 
for laser intensity fluctuations The detector was ac- 
curately calibrated and the software directly gives the 
number of atoms in the cold cloud. 



B. Continuous and pulsed regime 

When we turn on the ionic beam and send it on a hot 
neutralizer set at a constant temperature, we expect the 
spectroscopic signals to be well described by Eqs. © and 
© (Fig. [2]). The main advantage of this method is that 
the temperature of yttrium is kept constant during the 
measurement, and the vacuum conditions for the MOT 
are very stable. We analyzed Rb data taken in the con- 
tinuous regime for several temperatures. The comparison 
of the results for the fits with Eq. © and with Eq. © 
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showed an interesting feature, which allows us to esti- 
mate the impact of the implantation distribution on our 
results: while the two resulting curves are almost indis- 
tinguishable, assuming a sharp implantation distribution 
yields diffusion times that are systematically lower by a 
factor 0.7 with respect to the results obtained in the hy- 
pothesis of an extended initial distribution such as the 
one used by Melconian et at. This means that the ac- 
curacy for the measured diffusion coefficient is anyway 
limited at the 30% level by_the theoretical model, also 
for previous experiments 
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FIG. 2: Fluorescence signal from the Rb atoms desorbed by 
the neutralizer. The curve from Eq. @ (with an additional 
offset) was fitted to the experimental data. 

In order to enhance the signals, it is possible to operate 
in a pulsed mode. We let the ions accumulate in cold 
yttrium: at room temperature, the diffusion coefficient is 
very small, and the ions stay embedded in the neutralizer. 
When the neutralizer is turned on, the atoms are released 
all together in a time scale of the order of Td, according to 
Eq. ((T|) or In practice, for both Rb and Fr we usually 
kept the ionic current on for 600 s (which corresponds 
to more than three 210 Fr radioactive half-lives). This 
method gave excellent results to estimate the neutralizer 
release efficiency, and to determine the best operating 
temperature for our neutralizer. 



V. EXPERIMENTAL RESULTS 
A. Rubidium 

Most measurements have been done in the pulsed 
mode, because of the higher signal-to-noise ratio. We 
obtained a set of data for each neutralizer operating tem- 
perature: the main goal was to optimize the yttrium tem- 
perature in order to obtain a good release efficiency. 

The function G(G , S , r d , t ; t) = G + 5 /(l + (t - 



*o) / T d) 3 ^ 2 '', corresponding to the flux F Q (t) (Eq. was 
fitted to each acquired data set, with the free param- 
eters Go, So and Td (Fig. [3]). The offset Go is due to 
the presence of Rb impurities in yttrium: even with 
fresh neutralizers, never exposed to Rb current, we de- 
tect Rb vapor in the cell coming from the heated yttrium 
foil. Note that to is not left as a free parameter, be- 
cause it is totally correlated with the other parameters 
So and t^: for any value f , we can find the parameters 



S' = Sb/(1 + ^) 3 / 2 and r' d = r d (1 + ^), such 
that G(G ,So,T dl t ;t) = G(G 0) S ,T d ,t' ;t). Therefore, 
to is manually set to the time at which the signal begins 
rising, just after turning on the neutralizer. Since the 
neutralizer does not reach the operating temperature in- 
stantaneously, to is affected by an error corresponding to 
the signal rising time. 

The initial part of the curve critically depends on the 
way the neutralizer is turned on, and on the initial dis- 
tribution of ions inside yttrium (cf. Eqs. ([T]) and @). 
Therefore, we decided to consider data after a time i; n , 
corresponding to a signal which is equal to 3/4 of the 
maximum signal for the whole curve (Fig. [3]). 
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FIG. 3: Fluorescence signal from Rb neutralized ions in the 
pulsed mode and fit curve from Eq.Q, at a temperature of 
960 K. 

Release time data are consistent with Arrhenius' law 
for the diffusion coefficient D. Our experimental results 
are then fitted accordingly (Fig. [4]): 



Td = tiooo exp 



k n 1 

k B \T 1000K 



(7) 



where tiooo is the diffusion time at 1000 K; fit parameters 
are both the time tiooo and the activation energy E a . The 
choice of such a fit function lies in the fact that tiooo 
better describes the diffusion time in the temperature 
range of our measurements, while would be affected 
by a larger error, due to the extrapolation to infinite 
temperature. Results of fit are Tiooo = (6.1 ± 2.3) s and 
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FIG. 4: Diffusion time of Rb in yttrium as a function of tem- 
perature, fitted with the Arrhenius' function (Eq.0). 

E a = (1.8 ± 0.4) eV. We used the TRIM code to 
evaluate the mean implantation depth for Rb in Y at 
3 keV, which gives a value of 4.9 nm. From this data we 
can then deduce the diffusion coefficient at 1000 K: 

d 2 

£>iooo = ~ A = (1-0 ± 0.4) x 10~ 14 cmV 1 . 

4r iooo 

We also performed measurements on Rb using the MOT 
detection method, to prepare the Fr measurements and 
to test our diffusion model also for the MOT signal ev- 
lution. Experimental data and fit curve for Rb MOT 
in pulsed regime are shown in Fig. ([S]). A good agree- 
ment was found between data and model, and also with 
the results obtained from vapor fluorescence detection 
described above. The fit procedure used in this case is 
discussed in Section IVBl 




200 400 600 800 



Time (s) 

FIG. 5: Rubidium MOT in pulsed regime: experimental data 
and fit curve. Neutralizer temperature was about 960 K 



B. Francium 

The diffusion time of francium in yttrium has been 
measured using the magneto-optical trap. We operate in 
the pulsed mode in order to get higher signals: we use 
the same method as for Rb, with detection of trapped 
atoms. With a collisional rate C of about 0.5 s _1 , we 
checked that the condition /// <C C, W is verified for all 
the acquired data. 

As compared to vapor fluorescence data, MOT fluo- 
rescence is quite difficult to analyze in the pulsed mode, 
to deduce the diffusion time. When we turn on the neu- 
tralizer, francium atoms come out of yttrium along with 
many other impurities already contained in the neutral- 
izer (as Rb): vacuum suddenly deteriorates, and then 
slowly comes back to the equilibrium value. Since the co- 
efficient C in Eq. ([5]) strongly depends on vacuum quality, 
in the first part of the curve (t — to < Td) C depends on 
time and it is not easy to interpret the data. However, 
retaining only data with t — to 3> cannot be a solution, 
because in this regime So and Td are totally correlated (it 
is not possible to distinguish them in the fitting proce- 
dure) . Our compromise was to consider data after a time 
ii n that corresponds to a signal which is equal to half the 
maximum signal of the whole curve (Fig. [SJ. For these 
data vacuum effects are satisfactorily reduced and the 
order of magnitude is reasonable. 

It was not always easy to reach convergence in fits: in 
practice, we kept Sq and Td as free parameters, and we 
set the offset and to to a reasonable value for each fit. 
In some cases, Td and Sq were too correlated to obtain a 
good convergence: we had to repeat the fitting procedure 
with So set by hand to realistic values. The error for 
the obtained parameters includes several contributions: 
statistical error from the fit and uncertainties of Go and 
to- 

We saw in Section IV Al that the methodology was 
checked with Rb atoms. The agreement between MOT 
and vapor data demonstrates the validity of our measure- 
ments and of the fitting procedure. 

We expect to be able to improve our results. Because 
the presented measurements were initially intended to 
test the final release efficiency, in terms of number of 
atoms A*o rather than release time Td-, we preferred to 
begin with the pulsed method, which gives higher signals 
and higher signal-to-noise ratios. Nevertheless, in our 
limited beam-time we managed to acquire data in the 
continuous regime for one neutralizer temperature: we 
remind that in this regime we expect no vacuum-related 
problems. The analysis for the obtained curve gave very 
promising results, in agreement with the pulsed method. 

In Fig. [7] the results for the fitted parameter Td are 
reported as a function of temperature. We performed 
release time measurements using also the a decay rate 
signal, and we obtained similar results, at least as or- 
der of magnitude, thus confirming the consistency of our 
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FIG. 6: Trapped Fr fluorescence and fit curve according to 



Eq. © with T r = 3.633 x 1(T 
the CCD camera was set to 1 s. 



The exposure time of 
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FIG. 7: Diffusion time of Fr in Y as a function of temperature: 
experimental data and fit curve according to Arrhenius' law. 



method. Fit results according to Eq.([7|) are: 

r 1000 = (4.1 ± 2.5) s E a = (1.8 ± 1.1) eV 

An activation energy of E a = (1.59 ± 0.15) eV was 
measured for K in yttrium in Ref. [4|, so our results for 
both Fr and Rb are consistent with literature. We used 
the TRIM code 7] to evaluate the mean implantation 
depth for Fr in Y at 3 keV, which gives a value of 5.1 nm. 
The measured diffusion coefficient at 1000 K is then 



Diooo = (1.6 ±0.9) x 10 



" 14 cmV 1 



This is to our knowledge the first determination of diffu- 
sion time of francium ions in yttrium and the first time 
that diffusion parameters are derived from the study of 
the time evolution of the signal from a magneto-optical 
trap. 

According to Ref. @ it is possible to define a parame- 
ter e(T r ) = exp (— 2\/T r r^), which takes into account the 



loss of atoms due to nuclear decay during the diffusion 
towards surface. This is the so-called "release fraction" 
which gives the percentage of radioactive ions neutralized 
before the decay and then suitable for trapping. In other 
words, one could interpret s(T r ) as a sort of neutraliza- 
tion efficiency. From the diffusion time measurement we 
determine the release fraction as a function of tempera- 
ture as shown in Fig. [5] 
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FIG. 8: (Color online) Calculated release fraction vs tem- 
perature for isotope 210 (black curve) and 209 (blue curve), 
according to Arrhenius' law and E a and riooo given from our 
fit. Dashed curves give the uncertainty interval. 

Note that our measurements of the diffusion coeffi- 
cients for 210 Fr allow us to predict the release fraction 
for 209 Fr, which we also produce and trap, since we can 
reasonably suppose that diffusion features are identical 
for the two isotopes with negligible mass difference. Be- 
cause the apparatus was not optimized for the production 
of this isotope, we did not perform systematic measure- 
ments on 209 Fr. Nevertheless, our prediction of e provides 
an estimation of the optimal conditions for the 209 Fr neu- 
tralization: due to shorter lifetime of about 72 s, we ex- 
pect 209 Fr neutralization to be more efficient at higher 
temperatures with respect to 210 Fr (Fig. [5]). Analog con- 
siderations can be applied to 211 Fr: due to very similar 
radioactive lifetime, the curve calculated for 211 Fr is in- 
distinguishable from the 210 Fr one. 



VI. CONCLUSIONS 

We measured for the first time the diffusion param- 
eters of francium and rubidium ions in yttrium. The 
methodology that we developed takes advantage of laser 
spectroscopy techniques, and of the very high sensitivity 
offered by laser cooling and trapping. Since we directly 
measure diffusion times, our method can be considered 
complementary of release efficiency measurements by nu- 
clear decay. Another advantage is the possibility to mea- 
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sure diffusion coefficients aiso for stable elements. 
Experimental data are consistent with the diffusive model 
considered, demonstrating the validity and the potential 
useful application of our method. The diffusion coeffi- 
cient is determined together with an activation energy 



consistent with previous measurements. For Fr, a release 
fraction of more than 80% is obtained for the 210 isotope 
at 1050 K, while for the short-lived 209 isotope we predict 
a release fraction around 70% at the same temperature. 
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